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Abstract. In this paper we study a distance-dependent surface tension, defined as 
the free-energy cost to put metastable states at a given distance. This will be done 
in the framework of a disordered microscopic model with Kac interactions that can be 
solved in the mean- field limit. 



Surface Tension in Kac Glass Models 



2 



1. Introduction 

The introduction of point-to-set correlation functions [H El [3l HI |5] allowed important 
progresses in understanding the growth of static correlations in supercooled liquids near 
the glass transition. These non-standard correlation functions measure how deeply the 
effect of amorphous boundary conditions penetrates within a system. In order to intro- 
duce them, let us consider a large ensemble of interacting particles becoming glassy at 
low temperature. We assume that the liquid is trapped in a metastable state. We freeze 
the motion of all the particles outside a sphere of radius R. Then we let the particles 
inside the sphere free to move and eventually to rearrange in a different metastable 
state. The effect of the external particles is to create a pinning field favouring internal 
configurations which best match the frozen exterior. For small radius R, the effect of 
the pinning field on the interior of the sphere is strong. In that case the sphere remains 
in the same state. On the contrary, for large radius R, the effect becomes weak and 
the sphere can be found in a different state. Roughly speaking, a point-to-set correla- 
tion function measures the overlap between the initial state and the one reached after 
the rearrangement of the system. It has been found in numerical experiments that on 
lowering the temperature the effect of the amorphous boundary conditions propagates 
deeper into the region [3^, . 

Standard Random First Order Transition (RFOT) [6l [1] assumes that the 
competition between an entropy-rich state with high energy and an entropy-poor state 
with low energy, can explain the transition from high-overlap to low-overlap metastable 
states of the previous system, as the radius of the sphere is increased. As we are 
going to show, such a mechanism has to be reconsidered. In order to do this, let us 
consider, for simplicity, a Ising-like model described by an Hamiltonian H. We freeze a 
configuration S*" in a region A of the system. We study the thermodynamic considering 
only configurations S constrained to be close to S*" is A: 

Z[S"]=^e-^^[^]xA[5,^1 , (1) 
s 

where 

The thermodynamic average of an observable O of the system is obtained by averaging 

with constrained Boltzmann measure the configurations inside the sphere and with 

Boltzmann measure the configurations S": 

-[SH[s-] f,-ms] 
(O) = Y.^^J2xa[S,S-]——0{S) . (3) 

s ^ -I 

This average coincides with the usual thermodynamical one: -^Yls^~^^^^^^i^)- This 
simple fact has deep implications: in the case in which A is a sphere of radius R, on 
average, the energy per degree of freedom is independent of R. If, for typical choices of 
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the position of the sphere, one finds that two thermodynamic states coexist for a well 
defined value of R, they will have the same energy. Possible mechanisms for coexistence 
should therefore have a purely entropic origin [7]. 

In recent numerical experiments [8J the energy paid to put different metastable 
states in contact has been measured. The procedure is the following: freeze two states 
a and /3, exchange a sphere of the state a with a sphere of the state /3 and let the system 
evolve. Inspired by this idea, in the present work we want to introduce a different point- 
to-set correlation function defined as the free-energy cost to put different metastable 
states at distance /. In order to do that, we consider a sandwich-geometTj: two regions 
of the space divided by a box of width / and then freeze different metastable states at 
opposite sides of the box. Figure [H This system is well suited in order to be studied in 
the framework of a p-spin model with Kac interaction [5l |9l HDl [HI [121 113] . 

The paper is organized as follows: in Section [2] we introduce the model that we 
consider and the basic definitions; in Section [3] we briefiy illustrate how to obtain the 



free energy of the system; more details on these calculations can be found in Appendix A 



and in Appendix B; in Section H] we present our results and in Section [5] we draw our 



conclusions. 
2. The model 

We consider a finite-dimensional version of the spherical p-spin model, defined on a 
(i-dimensional cubic lattice A of linear size L, whose elementary degrees of freedom are 
spins G M with i G A. We introduce the interaction range 7^^ > and a non negative 
rapidly decreasing function iIj{x) normalized by: J d'^xtp{\x\) = 1. We define the local 
overlap of two configurations and 5*^ as: 

QsM^) = YY.'f'{l\^-j\)SjS] . (4) 

j&A 

We impose that configurations are subjected to the local spherical constraint: Q51 51 {i) = 
1 G A. We then introduce the finite-range p-spin Hamiltonian: 

Hp[S,J] = — Ji-i,...,ipSi^...Si^ (5) 
where the couplings Jii^,„^i are i.i.d. random variables with zero mean and variance: 



E 



, . . . , tp 

2 



J 

2l ,...,2p 



7 

fcGA 



is the interaction range since only variables located at vertices i and j such that 
< 7^ really interact. This also implies that the Hamiltonian is a random variable 
with zero mean and variance: 

E [H[S\J]H[S^J]] =Y^f{QsM^)) , (7) 
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where f{x) is a polynomial with positive coefficients, for example f{x) = x^, if we 
consider a pure p-spin model; in the following we will consider f{x) = j^x'^ + x^, where 
the quartic term assures to have a regular gradient expansion of the free-energy density. 

We analyze the model in the Kac limit: i^,7~^ — ?■ oo with L 3> 7~^, where the 
model can be solved by saddle-point approximation. 



a 
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a 


I 


f3 
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Figure 1. The sandwich-geometry for a system aa (left) and a/3 (right). 
B{1) is the central region, A~^{1) and A~{1) are the lateral ones. 



The box 



The sandwich-geometry is implemented by considering three regions of the lattice 
A: A'^{1), A^{1) and a box B{1), Figure [H In order to put the same or different states 
at opposites sides of the box, we introduce two different systems, that we call aa and 

af3: 

• system aa: we fix a configuration S*" drawn from the Boltzmann equilibrium 
measure. We consider the thermodynamic of configurations S constrained to be 
close to S*" both in A~^{1) and in A~{1); 

• system a/?: we fix two configurations S"" and drawn from the Boltzmann 
equilibrium measure. We consider the thermodynamic of configurations S 
constrained to be close to 5*" in A'^{1) and to S'^ in A~{1). 

We consider a system a/3. Let be O an observable of the system and q < 1. The 
constrained Boltzmann measure is: 

X n SiQs-sil) - g) n ^(Qs^si^ - q) (8) 

where j denotes integration over configurations satisfying the local spherical constraint. 
The partition function is: 

X n ^i.Qs-s{i) - g) n ^^Qsf^si^ - q) ■ (9) 

ieA- ieA+ 
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The symbol E represents the average over both the distribution of fixed configurations 
and S*^ and the disorder; the free energy of the system F^pil) is then: 



1e 

(3 



lnZ[^2+,^5_ 



(10) 



For a system aa, the constrained Boltzmann measure {■)aa{l) is obtained by imposing 
the constraint IljeA+uyi- ^{Qs'-s{i) - q); then Faa{l,T) = -^E [in Z[S';^+^^_]] 

As we will see in the following, F^^ll, T) and F^ail, T) can be calculated in the Kac 
limit, 7—7-0 taken after L — )■ oo. This allows us to measure the free-energy cost per 
unit area to put different metastable states at a distance /: 

Fal3{l,T) — Faa{l,T) 



Y{l,T) = lim lim 

7— >-0 L— >-oo 



11) 



this quantity can be interpreted as an effective, distance-dependent, surface tension. 



3. Calculations 



In the following we consider a system af3; a system aa can be treated in the same way. 
In order to calculate Fap, the average E can be taken by introducing replicas along the 
lines of [T2| 113] (more details on calculations can be found in Appendix B ). Integrals 
over spin variables are then treated for an (m + n) x (m + n) matrix order parameter 
<lab{i)- We rescale the position to define x = ij E [—L, L^, L = 7L to get: 

1 / r , -:KSa0{qab) 



^a/3(0 = lim 



[dqab]e~ 



(12) 



The dependency upon 7 is now completely explicit and, for 7 — )■ 0, the functional integral 
can be performed using the saddle-point method. We look for a replica symmetric saddle 



point 



ab \ 



This is characterized by three scalar functions Pi{x), P2{x) and q{x); pi 



and P2 are the local overlap between the constrained configuration and the reference 
configuration S°' and respectively and q is the local overlap of two constrained 



configurations when they belong to the same metastable state (see Appendix A for 
more details). Using this ansatz we obtain that Sapiqab) = ^ j J^apf^'^'x + 0{n'^), where: 

Ca^ix) = - ^[/(l) + 2/((^ * pi)(x)) + 2/((^ * p2){x)) - /((^ * q){x))] + 



1 

+ 2 



log(l - g(x)) 



pI{x) 



pI{x) 



q[x) 



1 ~ q{x) 



(13) 



with: 



{i^*q){x)= / d'^yipily - x\)q{y) 



(14) 



The constraint enforcing S to be close to S*" in A~{i) and to in A'^iJ) is fulfilled by 
setting pi{x) = g for x G A~(/) and P2{x) = g for x G A'^(i). We obtain Fa/sil) by evalu- 
ating the fields Pi{x), P2{x) and q{x) in the saddle point of the action 5°^ = / d'^xCapix). 
The resulting free energy will present an extensive part 0{L'^) which will be the same 
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for a system aa and for a system a/5. Then, in the calculation of the surface ten- 
sion T), the extensive part of the free energy will erase and contributions come 
only from the sub-leading order 0{L'^~^); the resulting form of the surface tension is 
F(/,T) = F„^(/,T) - F,„(/,T), where F„^(/,T) = ^ dxC^pix). 

We introduce a simplification in the Lagrangians: we expand the terms of the 
form /{{ip * in gradient of q{x) and we truncate to the second order obtaining 

— c/"(g(x))(Vg)^(x) where c = ^ J z'^ip{\z\)dz'^ (in our running example c = 1). 
We find the saddle-point fields iterating numerically the Euler-Lagrange equations of 

4. Results 

The system aa has been studied in spherical geometry we verified that in sandwich- 
geometry the behaviour does not change with respect to the spherical one. Two critical 
temperatures characterize the system: ~ 0.766287 and ~ 0.813526. 

Setting the temperature of the system T > T^, we find two lengths: lo{T) and 
$,diT), such that, for widths of the box 1 G [fo(T), ^^(T)], the action has two local 
minima. A minimum is characterized by a saddle-point field p{x) rapidly decaying 
to zero in the interior of the box; we name this low-overlap minimum. The other 
minimum is characterized by a saddle-point field p{x) everywhere large; we name this 
high-overlap minimum. For / > (/ < lo) only the low-(high-)overlap minimum exists. 
^s(T) is defined as the minimum value of I such that the low-overlap minimum is the 
global minimum of the action. The critical temperatures Tg and Td are defined as the 
temperature at which ^^(T) and C,d{T) respectively diverge. 

For a better comprehension, we present in Figure [3] the plot of the sub-extensive 
part of the free energy of high- (low-) overlap minimum F^^{i) {F^^{i)) divided by 
the size / for a system at a temperature < T < T^. C,s{T) is then the value 
of I where F^^il) and Fj^^il) cross. Then the global free energy of a system aa is 
F^4)=min[F^J),F^J)]. 

On the other hand, in the case of a system a/?, the action S^^ has always a single 
minimum. Profiles of the saddle-point field pi {x) can be seen in Figure [2l The sub- 
extensive part of the free energy of the unique minimum Fq,^(/)// for a temperature 
Ts < T < is also plotted in Figure [31 At all temperatures and values of / that we 
have studied, the sub-extensive part of the free energy of a system a/? Fapil) is close to 
the sub-extensive part of the low-overlap free energy of a system aa Faail), as can be 
seen in the inset of Figure [31 

In Figure[llwe follow the evolution of /-dependent surface tension Y{1, T) for systems 
at different temperatures T > Tg. We note that the static correlation length is{T) 
separates two regimes. For / < ^s(T), Y{l,T) has a power- law followed by a linear 



Surface Tension in Kac Glass Models 



7 



decrease. For / > is{T), as we see in the inset of Figure HI the decrease becomes 
exponential: 

r(/,T) ~ C e"'/', (15) 

with / weakly dependent on the temperature and showing no evident relation with 
^s- This shows that the surface tension Y(i, T) is sensibly different from zero only for 
^ ^ is- A similar result has been obtained in |T1]; in that case the interface free energy 
has been obtained changing the boundary conditions along one direction, from periodic 
to anti-periodic. 

Particular attention must be spent in the case T = T,. At T^, the static correlation 
length diverges. This means that the high-overlap minimum is the global minimum of 
the action iS°„ for all the values of /. We see in Figure H] that, for T approaching T^, the 
profile of F (/, T) takes the shape of a plateau. Consequently, at the critical temperature 
Ts, in the limit / — t- oo, the surface tension F (/, T^) does not fall to zero and takes a 
limiting value Y{Ts). Arguably, the value Y{T) is different from zero for temperatures 
T <Ts. 

According to phenomenological arguments jT], the static correlation length ^^(T) 
can be interpreted as the typical size of metastable states of a system at a temperature 
T. Following this idea, in a system afi we are freezing a patchwork of metastable states 
of size is{T) outside the box and letting the system free to rearrange inside the box. 
If the width of the box is larger than the typical metastable-state size, 1 ^ is{T)^ the 
system inside the box has enough space to rearrange in many different metastable states. 
On the contrary, when the width of the box is smaller than the metastable-state size, 
/ < since there is not enough space to create metastable states on the interior, the 
frozen states are in contact and then "repel" each other. This explains why the surface 
tension F (/, T) is significantly different from zero only for 1 < is(T) and why the overlap 
profiles pi{x) and P2{x) between frozen metastable states and the interior of the box 
decrease faster for small boxes. At the critical temperature Tg the size of metastable 
states diverges. Consequently, the surface tension takes a finite value also in the limit 

/ — 7- OO. 



Other observables of the system have been considered. We studied the internal 
energy U. We verified that for a system aa the high-overlap and the low-overlap phases 
have the same energy, as motivated in Section [H In Figure [5] we follow the evolution of 
Uajiil) — Uaa{l) foT different temperatures of the system. A detailed derivation of this 



quantity can be found in Appendix B In this case, we note a power-law followed by a 
an exponential decrease. 

We also computed the configurational entropy S as a function of the size / of the 
box. Figure [6l For a system aa only the low-overlap phase presents a configurational 
entropy Ti^a different from zero. As noticed in [5j, for I < 1^^^^ the replica-symmetric 
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solution is incorrect since it gives a negative entropy. We found that the same is true for 
a system a/3. In the inset of Figure [6] we plot the difference between the configurational 
entropy of the two systems. We note that this quantity is a decreasing function of 
the size / of the system. This is consistent with the observation that the system loses 
memory of the frozen exterior for large sizes of the box. 



5. Conclusions 



In this paper we have studied a distance-dependent surface tension, defined as the 
free-energy cost to put metastable states at a given distance. This has been done in the 
framework of a disordered microscopic model with Kac interactions that can be solved in 
the mean-field limit. We have found that the surface tension is sensibly different from 
zero only for distances between metastable states smaller than the static correlation 
length of the system. A description of this behaviour in terms of a phenomenological 
droplet argument has been proposed. Other observables, like the internal energy and the 
configurational entropy, has been studied. The behaviour of the configurational entropy 
allowed to identify under which size the replica-symmetric ansatz becomes incorrect and 
a 1-RSB solution must be considered. 
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Appendix A. 



We give an explicit formulation of the overlap matrix in the replica symmetric ansatz 
g^jf (x); the overlap matrix is an {m + n) x [m + n) matrix with m and n real numbers; 
if we take r and n integer, we can visualize the matrix in the following way: 

A B 



RS 



(A.i; 



The n X n matrix C is the overlap matrix between configuration that are taken with 
constrained Boltzmann measure and subjected to local spherical constraint; the replica 
symmetric ansatz imposes that Cab = for all a 7^ 6 and the spherical constraint 
that Caa = 1; then C can be written in the form: 

1 q{x) q{x) 



C 



q{x) 1 q{x) 



q[x) 
q{x) 



(A.2) 



q{x) q{x) q{x) ... 1 

The m X m A matrix is the overlap matrix between configuration that are taken with 
Boltzmann measure and subjected to local spherical constraint; we impose the out of 
diagonal elements equal to zero, then we obtain that A is the identity matrix: A = 1. 
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The m X n B matrix is the overlap matrix between configuration that are taken 
with Boltzmann measure and configuration that are taken with constrained Bohzmann 
measure; we impose ah the elements of this matrix equal to zero, except the last two 
lines that are equal to pi{x) and to P2{x); then A can be written in the form: 



B 






pi{x) pi{x) 

P2{x) P2{x) 






pi{x) 

P2{x) 



(A.3) 



Appendix B. 

The internal energy of a system a/? is: 

C/c,^(0=E[(//[s,J]),^(0]=E 



d_ 

"dp 



We give an explicit derivation of Uaj3{l)', similar calculations allow to obtain the free 
energy Fap{l). We introduce two different temperatures and ^2 and n and m replicas 
of the system; 



Uo.^{l) = E 



dfi 



E 



dp2 



n,n->-0 n \ dp2 



E 



Z™-2[/3i] j ds'^ds^ exp [-^(^[s", J] + H[s^, J])] Z''[s'', 3^, 82 



where: 



m+n 

n 

a=m+l 



ieA+ ieA- 

Then performing the expectation value over the disorder, the derivative and reimposing 
equal the temperatures we obtain: 



„ ,. m+n 

lim ^ / TT ds"Cexp 

a=l 



ieA l<a,6<?i 



E /(Qaf(O) + ^ E fiQabii)) + J E fiQabii)) 



a,b£C 



a,bGB 



Integrals over the spin variables are then traded for an (m + n) x (m + n) matrix order 
parameter Qabii)- 



Ua^= lim — / IT TT 'lab{i)Cexp 

v^-^A^I 1 



ieA a,b=l 



ElC E /(MO)) + ^logdetQ(i) 



l<a,b<r 



E 

ieA 



E f(9ab{i)) + l E + J E /(^'^''W) 



a.fceC 



o.fceB 
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a,beC 



a,beB.BT 



X exp 



4l / ( Y E /((^ * + i log det q{x) 



c . 



Using the rephca symmetric matrix presented in Appendix A we obtain: 



Uaf3 = lim — / [dg][dpi][dp2] 



n / dx'H{x) + o{n ) 





n f 


exp 









where: 



■H{x) = 1 + /((V + /((V *P2)(x)) - /((V * q){x)) ; 



^(^) = ? [1 + 2/((V^ *pi)(x)) + 2f{{i^^p2){x)) - /((V. * + 



log(l - q{x)) 



pI{x) +pI{x) -q{x) 



1 - q{x) 



We evaluate the action = J d'^xCa/3{x) in the saddle point fields pi, p2 and q and 
we obtain that: 



f/„^(/) = -/3 / dxn{x) . 



(B.l) 
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Figure 3. Main figure: Plot of tlie sub-extensive part of tlie free energy divided by 
the size as a function of / for a system at a temperature T = 0.7874 of: high-overlap 
minimum of a system aa, F^^{l)/l; low-overlap minimum of a system aa, F^^(l)/l; 
unique minimum of a system a/3, Fai3(l)/l- The static correlation length is pointed 
out. Using this scale F^^{1) and Foc0{l) are indistinguishable. Inset: the difference 
F^^(l) — Fap(l) in logarithmic scale. 
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Figure 4. Plot of Y{1, T) for different temperatures as a function of the width of the 
box i. We remember that « 0.766287 and « 0.813526. 
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Figure 5. Plot of Uap(l) — Uaa{l) for different temperatures as a function of the width 
of the box /. We remember that w 0.766287 and f« 0.813526. 
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